Optimal control of circuit quantum electrodynamics in one and two dimensions 
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Optimal control can be used to significantly improve multi-qubit gates in quantum information 
processing hardware architectures based on superconducting circuit quantum electrodynamics. We 
apply this approach not only to dispersive gates of two qubits inside a cavity, but, more generally, 
to architectures based on two-dimensional arrays of cavities and qubits. For high-fidelity gate 
operations, simultaneous evolutions of controls and couplings in the two coupling dimensions of 
cavity grids are shown to be significantly faster than conventional sequential implementations. Even 
under experimentally realistic conditions speedups by a factor of three can be gained. The methods 
immediately scale to large grids and indirect gates between arbitrary pairs of qubits on the grid. 
They are anticipated to be paradigmatic for 2D arrays and lattices of controllable qubits. 

PACS numbers: 03.67.Lx, 85.25.-j, 82.56. Jn 



I. INTRODUCTION 

Progress towards the goal of scalable quantum infor- 
mation processing is currently concentrated in physical 
systems that live at the intersection between quantum 
optics and solid state physics. One of the most promis- 
ing contenders are superconducting circuits that couple 
qubits and microwave resonators, a new field now known 
as circuit quantum electrodynamics (QED). In the de- 
velopment of this field, early suggestions to implement 
the Jaynes-Cummings model in the solid stated - — were 
followed by a seminal proposal^ to employ on-chip mi- 
crowave resonators and couple them to artifical atoms 
in the form of superconducting qubits. Experimental 
realizations soon followed^, creating a solid-state ana- 
logue of conventional optical cavity QED&. The tight 
confinement of the field mode and the large electric 
dipolc moment of the 'atom' yield extraordinary cou- 
pling strengths, which has led to a variety of experi- 
mental achievements, including: The Jaynes-Cummings 
model in the strong-coupling regime^i^, Rabi and Ram- 
sey oscillations and dispersive qubit readout*2ii£, gener- 
ation of single photons^ and Fock state a 12 ' 13 , cavity- 
mediated coupling of two qubit o 14 ' 15 , setups with three 
qubits^, Berry's phased, and the measurement of the 
photon number distribution^. 

The recent experimental progress in creating mi- 
crowave circuits with multiple qubits coupled to res- 
onators establishes the need for efficient, high-fidelity 
multi-qubit quantum gates and motivates the search for 
advanced architectures for quantum information process- 
ing on the chip. The most elementary situation to con- 
sider is two qubits inside a cavity (transmission line res- 
onator). The mere presence of the cavity induces a flip- 
flop (XY) type interaction between the qubits, which 
may and has been used for entangling gate operations^. 
The XY interaction directly produces an iswAP two- 



qubit gate. Other gates (like the CNOt) have to be syn- 
thesized. While it is well-known how to do so using a 
sequence of iswAP and single-qubit gate o 19 i 20 , there is a 
lot of room for improvement in constructing faster gates 
even for this basic situation. One way to go is to use 
resonant two-qubit gate a 21 ' 22 . The other approach is to 
keep the robust dispersive XY interaction and to explore 
better pulse sequences using the tools of optimal control 
theory. This is the approach we will follow in this paper. 

When going beyond two qubits and connecting many 
qubits into a quantum processor on a chip, it is crucial to 
abandon linear arrays and to extend the setup into the 
second dimension. While there have been a number of 
schemes for doing so in nearest-neighbor coupled 2D ar- 
rays, the presence of global coupling between qubits via 
resonators adds a new feature that has to be explored. 
Perhaps the most straightforward route, recently intro- 
duced by some of ns^, is to create a two-dimensional 
'cavity grid' of transmission line resonators arranged in 
columns and rows (Fig.[T|). The qubits are located at the 
intersections, such that each qubit feels the microwave 
field of two cavities. In this way, qubits can be addressed 
and coupled. An interesting feature specific to such a 
global coupling architecture is the fact that two qubits 
placed anywhere on the grid can be coupled via a third 
qubit sitting at the intersection of two resonators, with 
an overhead that does not grow with system size. Again, 
while a sequential protocol for this 'coupling around the 
corner' is known 2 ^, we may ask for possible speedups ac- 
quired by more sophisticated pulse schemes. 

Principles of optimal quantum control^ 3 , are currently 
establishing themselves as indispensible tools to steer 
quantum systems in a robust, time-optimal or loss- 
avoiding way 2 ^. Earlier examples of numerical optimal 
control can be found in Ref. [25|. Based on the gradient 
ascent approach grapes 2 ^, our toolbox of optimal control 
techniques now allows the synthesis of quantum gates 
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Figure 1: (Color online) The superconducting cavity grid^S, 
with two layers of vertical (bottom) and horizontal (top layer) 
transmission line resonators, coupled to qubits (small red 
squares). Two-qubit gates between qubits 1 and 3 are me- 
diated indirectly via qubit 2, employing the dispersive inter- 
action inside the two highlighted resonators. 



in a time-optimal2l or relaxation-optimized way, where 
the open system may come in a Markovian^ or a non- 
Markovian setting^. 

In particular in superconducting qubits (recently re- 
viewed in Ref. [30h . the issue of improving qubit gates 
by optimal control techniques is currently attracting in- 
creased attention (see Refs. l3lM36l for some of the ap- 
proaches). Here, we take these ideas further by opti- 
mizing multi-qubit gates in the setting of circuit quan- 
tum electrodynamics. We will consider both the stan- 
dard two-qubit setup, as well as the three-qubit config- 
uration where interactions exist only between qubits 1 
and 2 and qubits 2 and 3, while the gate is to be applied 
between qubits 1 and 3. Note that the latter situation is 
of far more general relevance than the cavity grid only. 
In fact, it will become important whenever full global 
coupling (all qubits couple to all others) is not available, 
which is the expected situation for multi-qubit circuits. 
Optimal control techniques have already been success- 
fully applied to other three-qubit architectures, in both 
analytical^— , and numerical^ approaches. 

The remainder of the paper is organized as follows: We 
first provide more background on the cavity grid archi- 
tecture and our optimal control methods in Section [TXJ 
Afterwards, we go through a sequence of models of in- 
creasing sophistication, for which we discuss the results 
found using optimal control. We look at two qubits in 
one cavity and three qubits in two cavities (the cavity 
grid situation). At first, in Section [ill! wc consider both 
of these situations in the setting of an 'idealized model' 
where we assume local control over all the qubits. Later 
on in Section IIV1 we introduce a 'realistic model', in 
which we take into account that a microwave pulse ap- 
plied through a resonator will couple to both qubits, such 
that some amount of local control is lost. In addition, 
in the 'realistic model', we will consider restrictions on 



the pulse amplitudes. In all cases, we are interested in 
finding the minimum time at which the optimized pulse 
sequence has fidelity unity, and we extract the speedup 
vs. the known, sequential approach. We will discuss 
typical control sequences and also present the evolution 
of entanglement between the qubits during an optimized 
three-qubit scheme. 

II. METHODS & SETUP 
A. The cavity grid 

We briefly recount the basic setup of the two- 
dimensional cavity grid2£. The idea is to have a grid 
of transmission line resonators (horizontal and vertical, 
in two different layers), and to place qubits at the inter- 
sections, as depicted in Fig. [T] Each qubit thus feels the 
microwave field of the two cavities crossing at its loca- 
tion. It thus couples directly, via the standard dispersive 
XY-type interaction erf erf J ra\a\ (see Section ITO]) , to all 
of the qubits in both of these cavities. Qubit frequencies 
are chosen distinct to provide for individual addressabil- 
ity. iswAP gates can then be implemented by bringing 
two qubits into mutual resonance (still detuned from the 
coupling cavity), evolving for an appropriate time, and 
bringing them out of resonance again. For an n x n ar- 
ray of n 2 qubits, only n different frequencies are needed, 
which drastically eases the restrictions for the size of the 
array vs. the one-dimensional case. 

A crucial question, however, is how to couple two 
qubits that are not part of the same cavity. In nearest- 
neighbor coupled arrays, this would require a number of 
intermediate SWAP operations that grows with the size 
of the array. In contrast, the overhead remains constant 
for the case of the cavity grid. Assuming direct couplings 
exist between 1 — 2 and 2 — 3, then any two-qubit gate 
between 1 and 3 can be implemented by first swapping 
the quantum information from 1 to 2, then performing 
the desired operation, and finally swapping back. Thus, 
since any two qubits will share a third qubit at the cor- 
ner of the cavities to which they couple, this provides a 
general, distance-independent coupling scheme. 

However, one of the disadvantages of such a straight- 
forward pulse sequence is that the SWAP gate itself is not 
elementary with respect to the XY interaction, which 
only yields an iswAP. Unfortunately, three iSWAPs are 
needed to compose one SWAP, of which two are necessary 
for the 'coupling around the corner' approach. Thus, this 
very important operation in the cavity grid architecture 
lends itself naturally to potentially significant improve- 
ments via optimal control techniques. 

B. Our optimal control approach 

The general framework in optimal control is to max- 
imise a figure of merit subject to steering a dynamic sys- 
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tern according to its equation of motion under experi- 
mentally admissible controls. In quantum information 
processing, a convenient figure of merit often chosen is 
the trace fidelity of the gate U(T) actually synthesised 
under the available controls at some final time T with 
respect to the desired target unitary gate [/target- For n 
qubits and setting N := 2™ the (squared) trace fidelity is 

f* := |i tr{uL got U(T)}\ 2 . 

Now, for closed quantum systems Schrodinger's equation 
of motion lifted to operator form reads 

u = -mu, 

where the total Hamiltonian H is composed of the non- 
switchable drift term and control terms Hj governed 
by piecewise constant control amplitudes Uj(tk) for tk £ 
[0, T] so as to give 

H(t k ) ^Hd + YjUjit^Hj. 

3 

Dividing the total time T into M intervals At k (of piece- 
wise constant controls) with T = J^k Aifc, the gate 
U(T) = Um Um-i ■ ■ ■ Uk ■ ■ ■ U% Ui is made of components 
Uk = cxp{— iH(tk)Atk}- Then the derivatives 

Rc ^ tr{[/t argot U M U M -i ■ • ■ U k+1 x 

(iAttHjUjUk-t ■ ■ ■ U 2 U t } ti{U^ get U(T)}*) 

can (with appropriate step size a r > 0) readliy be used 
for recursive gradient schemes like 

representing the simple setting of steepest ascent. Like- 
wise, conjugate gradient or Newton methods can be 
implemente d 26 ' 41 . Here we used the LBFGS variant of 
a quasi-Newton approach, as sketched in Appendix [A] 
We refer to algorithms of this kind as Gradient Ascent 
Pulse Engineering^ 6 - (grape) algorithms. 

Spin and pseudo-spin systems are a particularly pow- 
erful paradigm of quantum systems, in particular when 
they are fully operator controllable, i.e. universal. For 
this to be the case the drift and the control Hamilto- 
nians have to generate the full n-qubit unitary algebra 
su(2") by way of commutation 4 ^—. A simple example 
exploited in the STIRAP 4 ^ scheme is the Lie algebra of 
local control su(2) being generated by the Pauli matrices 
a x (pulse) and a z (detuning), whose commutator gives a y 
and thus introduces phase sensitivity. In the instance of 
the (realistic) two-qubit Hamiltonian examined in Sec- 
tion [TV] (Eqn. [3]), only different phase shifts Ai ^ A 2 
ensure full controllability. In practice for universality it 
suffices that (i) all qubits can be addressed selectively 
and (ii) that they form an arbitrary connected graph of 
Ising-typc coupling interactions. More recent analyses 



revealed that even partially collective controls maintain 
universality in different types of Ising or Heisenberg cou- 
pled systems as long as the hardware architecture gives 
rise to Hamiltonians with no symmetries^ 6 -. 

III. IDEALIZED MODEL 

In order to establish the lower limits on gate times 
in this scheme, let us first consider an idealized model 
where the control fields are unrestricted. A model which 
respects more closely the limitations in current experi- 
ments will be considered in the following section. After 
adiabatic elimination of the cavity mode 4 -, the effective 
qubit-qubit interaction Hamiltonian is of the form 

Hmt = IX J (crt<72 + °T CT 2~) 
7rJ 

= ~2 ( a ? a 2 + °1°2) J 

where J is an effective coupling constant determined by 
the qubit-cavity couplings and detunings. Evolution un- 
der iJj nt for a time of T = -j yields the so-called iswAP 
operationi^: 

/ 1 \ 

exp { _i ^ 5nt } = -i ~o ' (2) 

V 1 / 

a universal two-qubit gate which can be considered the 
'natural' gate of the coupling interaction. 

A. Two qubits in a cavity 

If the two coupled qubits are individually addressable 
by resonant microwave fields of tunable amplitude and 
phase, the total Hamiltonian in a frame rotating with 
the driving fields is 

2 

HgL® = tf int +X>(Q?(t)af + Q?(tK), 

i=l 

under the assumption that the two qubits are in reso- 
nance (i.e. they are set to the same frequency, but dis- 
tinct from the cavity frequency). Note that our simula- 
tions work within the rotating-wavc approximation and 
assume a sufficiently detuned cavity that has already 
been eliminated. Thus, we neglect both the Bloch-Siegcrt 
shift that would arise for extremely strong driving, as well 
as any AC Stark shift due to a strong cavity population. 
In order to make the bilinear control form of the Hamil- 
tonian more explicit, the microwave fields are specified 
in terms of real and imaginary parts fif (t) and Of (t) , re- 
spectively, rather than amplitude and phase. Through- 
out this article the controls and J-coupling are normal- 
ized as frequencies rather than angular frequencies, with 



4 



the 2-7T factors written explicitly in the Hamiltonians (and 
with h = 1). 

One approach to implement a general two-qubit gate is 
to decompose it into a sequence of iswAP gates and local 
operations, as discussed in Rcfs. IT9II20I . For example, a 
CNOT gate can be created from two iSWAPs, while a SWAP 
gate requires three. We refer to this as the 'sequential' 
approach. In this section we assume that local operations 
can be performed in a negligible time compared to the 
time required by the coupling evolution. Time-optimal 
pulse sequences for an arbitrary two-qubit gate can then 
be determined analytically via the Cartan decomposition 
of SU(4) as in Refs. |47||48| . The iswAP implementation 
suggested in Eqn. ([2]) is, unsurprisingly, already time- 
optimal. Time-optimal pulse sequences for the SWAP and 
CNOT arc provided in Fig. [2] A comparison of the times 
required by the different schemes is given in Table Q] - we 
find that even in this simple case the swap and CNOT 
can be sped up by a factor of 2. 
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Figure 2: Analytical pulse sequences for time-optimal imple- 
mentations of two-qubit gates: (a) the CNOT gate, where 1 is 
the control qubit and 2 is the target qubit, and (b) the SWAP 
gate. 



B. Three qubits in two cavities 

We now consider two qubits, each in a separate cavity, 
which arc coupled indirectly via an additional 'mediator' 
qubit placed at the intersection of the cavities. If local 
controls on all three qubits are available, the Hamiltonian 



(3) 
deal 



ttJ 

— (of of + a\a\ + of of + <%a\) 



+ 5>(n?(t)of + n?(tK). 



Gates can be implemented between the indirectly cou- 
pled qubits (1 and 3) in the sequential scheme via the 
SWAP operation, as depicted in Fig. [3] For example, an 
iSWAP between qubits 1 and 3 could be implemented as 
a sequence of seven two-qubit iSWAPs. 

However, these indirect two-qubit gates embedded in 
a three-qubit system can be implemented considerably 
faster using optimized controls. The analytical methods 
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arbitrary two-qubit gate 

Figure 3: The standard decomposition of an indirect two- 
qubit gate into direct two-qubit gates via the SWAP operation. 



used for determining time-optimal two-qubit gates can- 
not be applied here; instead we use the numerical tech- 
niques outlined in Section [XXJ For a fixed gate time, an 
initial pulse is chosen and iterated until the GRAPE al- 
gorithm converges to a maximum of the fidelity This 
procedure is repeated for a range of different gate times, 
allowing us to estimate the minimal time. Further de- 
tails about the numerics are included in Appendix [A] A 
plot of maximum fidelity vs. gate time for the case of 
an isWAPi3 gate is shown in Fig. |U In this case we find 
that a time of 1/ J is required to reach the threshold fi- 
delity. Minimal times for other indirect two-qubit gates 
are similarly calculated and the results are included in 
Table HI alongside the times required by the correspond- 
ing sequential schemes of decomposition into two-qubit 
iSWAPs. 





Figure 4: (a) Maximum achievable fidelity as a function 
of pulse duration in the three-qubit idealized model for an 
iSWAPi3 gate, (b) On a logarithmic scale we observe a sharp 
convergence to the threshold fidelity of 1 — 10 -5 , where the 
algorithm terminates. 

In order to illustrate how the optimized indirect 
two-qubit gates differ from the sequential schemes, we 
can consider the entanglement between directly coupled 
qubits over the time interval during which the controls 
are applied. For this we use the logarithmic negativity^, 
defined as 

E N {p)=\o^\\p TA \\ x 

where r a is the partial transpose, is the trace norm, 
and p is the reduced density matrix of the two-qubit sub- 
system. We choose the initial state 1 100) and apply the 
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Gate 


T soq (1/J) 


TV (1/J) 


speedup factor 


iSWAPi2 


0.5 


0.5 


_ 


CNOT12 


1.0 


0.5 


2 


SWAP12 


1.5 


0.75 


2 


iSWAPi 3 


3.5 


1.00* 


3.50 


CNOT13 


2.0 


1.00* 


2.00 


SWAP13 


4.5 


1.15* 


3.91 



Table I: Implementation times for a selection of direct and 
indirect two-qubit gates in the idealized model: T seq is the 
time required by decomposing the gate into two-qubit iSWAPs; 
Topt is the time required by the optimal control sequence. The 
times marked with an asterisk are determined numerically as 
the shortest times in which the GRAPE algorithm can reach a 
fidelity of 1 — 10 , with time resolution 0.05/ J. The time 
of 2.0/ J for the sequential implementation of a CNOT13 is a 
special case, where the two swaps in Fig. [3] can be replaced 
by iSWAPs^. 



iSWAPi3 operation while observing the entanglement be- 
tween qubit pairs 1-2 and 2-3, illustrated in Fig. [5] In 
the sequential scheme the mediator qubit is entangled ei- 
ther with qubit 1 or qubit 3. In the optimized case, as 
one might expect, the mediator qubit is simultaneously 
entangled with both. 




t (1/J) t (1/J) 

(a) (b) 

Figure 5: (Color online) Logarithmic negativity between qubit 
pairs 1-2 (blue line) and 2-3 (red, dashed line) for (a) sequen- 
tial and (b) optimized implementations of an iSWAPi3. 



IV. REALISTIC MODEL 

Allowing for unrestricted x and y control on each qubit 
yields lower bounds on what implementation times are 
possible. However we would also like to consider a re- 
stricted model of coupled superconducting qubits which 
is more feasible in current experiments. We allow for in- 
dividual tuning of the qubit resonance frequencies (z con- 
trols), but restrict ourselves to a single microwave field 
(x control) per cavity, where the microwave field is no 
longer required to have tunable phase. 



A. Two qubits with restricted controls 

Under these structural restrictions, the corresponding 
two-qubit Hamiltonian is 

2 

#reilW = ffint + E^K+^WK + ^f) (3) 

i=l 

where Q(t) is the amplitude of the microwave field and 
Aj (t) are the detunings of the qubit frequencies from the 
microwave carrier frequency. We consider two possible 
cases for the control functions: (i) the controls are unre- 
stricted, or (ii) the controls are restricted to the following 
ranges: 

|Ai(t)| < A max = 1000MHz 
|fi(t)| < a„ ax = 50MHz, (4) 

where the coupling constant was taken to be J = 21MHz 
(as in Ref. I20I ) in our numerical examples. Furthermore, 
in case (ii) we require that the controls start and end at 
zero with a maximum rise-time of 4ns, which should be 
feasible in current experiments 2 ^. In case (i) the results 
from Section IIII Al still apply - we need only to rewrite 
the local x and y pulses in terms of our new controls. 
For instance a 90° x-rotation on the first qubit can be 
decomposed as 

i?*(90°) = fl|(-180°) i??. 2 (45°) R z 2 (180°) i^ >2 (45°) 

and the other local x and y pulses can be similarly de- 
composed. Thus, the two-qubit times in Table|T]also hold 
for this case. In particular, this means that case (i) also 
yields the optimum times for an experimental setup with 
both x and y pulses available (variable phase of the mi- 
crowave drive) and without restrictions on the local gate 
controls. In case (ii ) the analytical methods are no longer 
applicable, as they require that local rotations can be ap- 
plied in negligible time. The isWAP can of course still be 
implemented by simply evolving under the coupling, but 
to find time-optimal implementations for other two-qubit 
gates we again apply the GRAPE algorithm. This can be 
forced to optimize only over controls in a certain ampli- 
tude or bandwidth rang o 26 i 29 . Fig.[6]contains the fidelity 
vs. pulse duration curves for two-qubit SWAP and CNOT 
gates. Examples of the optimized controls obtained by 
the GRAPE algorithm for the CNOT gate are provided in 
Fig. [7] 

Observe that the CNOT gate is self-inverse and the two- 
qubit Hamiltonian in Eqn. is real and symmetric. 
As we described earlier—, in such systems there may 
be palindromic control sequences as in Fig. [7^,. Their 
practical advantage lies in the fact that they may be 
synthesized on LCR terminals only using capacitances 
(C) and inductances (L) and no resistive elements (R) 
thus avoiding losses. In contrast, since the iSWAP is only 
a fourth root of the identity, it is no longer self-inverse, 
and therefore palindromic controls are not to be expected 
(compare, e.g., Fig. El) 
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(a) (b) 

Figure 6: (Color online) Maximum fidelity as a function of 
pulse duration in the two-qubit realistic model for (a) a CNOT 
gate, and (b) a SWAP gate. Maxima obtained with no restric- 
tions on the controls are shown in blue (•), while those ob- 
tained under the restrictions in Eqn. Q are shown in red (o). 
The arrows indicate the minimal times at which the threshold 
fidelity is achieved. 




0.1 0.2 0.3 0.4 0.5 0.2 0.4 0.6 0.8 



t (1/J) t (1/J) 

(a) (b) 

Figure 7: (Color online) Sample controls obtained by the 
GRAPE algorithm for the minimal-time implementation of a 
CNOT in the two-qubit realistic model: Ai (red, dotted line), 
A2 (green, dashed line), and Q (blue, solid line) with (a) un- 
restricted controls, and (b) restricted controls. Note that the 
time course of the controls in (a) is perfectly palindromic as 
the controls may equally well be read forward and backward 
in time. 



A comparison of times in the sequential and optimized 
schemes under the restrictions in Eqn. (0]) is provided in 
Table [TTJ The times in the sequential scheme have in- 
creased, as each local 90° x-rotation now requires a time 
of 0.25/fi max = 0.105/ J. The times required by the op- 
timized schemes also increase, but substantial spcedups 
are still possible. 




0.8 1 1.2 1.4 1.6 1.8 0.8 1 1.2 1.4 1.6 1.8 



T (1/J) T (1/J) 

(a) (b) 

Figure 8: (Color online) Maximum fidelity as a function of 
pulse duration in the three-qubit realistic model for (a) a 
CNOT13 gate, and (b) an iSWAPi3 gate. Maxima obtained 
with no restrictions on the controls are shown in blue (•), 
while those obtained under the restrictions in Eqn. Q are 
shown in red (o). The arrows indicate the minimal times at 
which the threshold fidelity is achieved. 
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B. Three qubits with restricted controls 

For three qubits coupled via two cavities we allow for 
three local z controls and two x controls, with the Hamil- 
tonian 

+ 7rfi£,(t) (of + of) + 7Tfif 3 (t) (of + of) 

3 

+ 5>0?(i)of. 

i=l 

Again we determine optimized controls numerically; fi- 
delity vs. pulse duration curves are shown in Fig.[8l while 
sample optimized controls for the restricted case (ii) are 
shown in Fig. [5J Only the x restriction plays a role 
here as the z restriction is an order of magnitude larger. 



Figure 9: (Color online) Sample controls to implement an 
iSWAPi3 gate in the three-qubit realistic model with restric- 
tions (ii) in place: (a) Qf 2 (green, dashed line), O23 (blue, 
solid line), (b) Of (blue, solid line), Of (green, dashed line), 
03 (red, dotted line). 



V. CONCLUSIONS 

We have demonstrated how optimal control methods 
provide fast high-fidelity quantum gates for coupled su- 
perconducting qubits. In contrast to conventional ap- 
proaches that make use of the coupling evolutions se- 
quentially (i.e., along one dimension at a time), numer- 
ical optimal control exploits the coupling dimensions si- 
multaneously thereby gaining significant spcedups. In 
particular, our numerical method provides constructive 
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Gate 


T soq (1/J) 


TV (1/J) 


speedup factor 


iSWAPi2 


0.50 


0.50 


_ 


CNOT12 


1.21 


0.90 


1.34 


SWAP12 


1.82 


0.80 


2.28 


iSWAPi 3 


4.13 


1.40 


2.95 


CNOT13 


2.21 


1.40 


1.58 



Table II: Implementation times for a selection of direct and 
indirect two-qubit gates in the realistic model with the con- 
trol amplitudes restricted as described in case (ii). T aeq is the 
time required by decomposing the gate into two-qubit iSWAPs 
and local operations; T op t is the time required by the numer- 
ically optimized pulse to reach a fidelity of 1 — 10~ 3 . The 
particular values for the minimal times simply result from 
our choice for the maximum amplitude relative to J. 

controls under realistic experimental conditions, such as 

(i) power and rise-time limits in the control amplitudes, 

(ii) additional individual dctunings on each qubit, and 
(Hi) lack of phase switching and restriction to controls 
affecting qubits jointly. We showed how the latter two 
complement one another. Finally, our approach to pro- 
viding optimized controls tailored to the hardware of 2D 
cavity grids scales to (2D) arrays of many qubits and ar- 
bitrary gate operations between two generic qubits picked 
from such grids. It is therefore anticipated to find wide 
application in similar architectures of large-scale quan- 
tum systems used as quantum simulators, processors, or 
storage devices. 
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Appendix A: Numerical details 

The gradients introduced in Section III Bl form the ba- 
sis of our numerical optimization algorithm, but there is 
still some freedom in how they can be used to update the 
controls. The simplest approach of adding the gradients 
directly to the controls with some positive stepsize, as in 
Eqn. ([T]), is typically not the most numerically efficient 
one. Here we replace the conjugate-gradient updates of 
our standard routine^ by a Newton method which makes 
use of the second derivative of the quality function (the 
Hessian matrix). Because we chose to divide our evolu- 
tion into M = 256 constant intervals, this results in a 
256 x 256 Hessian matrix, direct computation of which 
would be inefficient. Instead we use the standard limited- 
memory Broyden-Fletcher-Goldfarb-Shanno (lbfgs) al- 
gorithm to approximate the Hessian matrix by analyzing 
successive gradient vector a 41 ' 50 . 

Inherent to these gradient-based methods is that they 
result only in a local maximum of the fidelity. To increase 
confidence that optima close to the global maximum will 
be obtained, we can sample from a range of random initial 
controls. For the results presented in this paper we used 
the following procedure: 

(i) Randomly generate 50 initial controls and optimize 
each in 100 iterations of the GRAPE algorithm. 

(ii) Select the highest 10 fidelities and iterate each a 
further 500 times. 

(ra) Select the highest 2 fidelities and iterate each a fur- 
ther 1000 times. 

This procedure is applied for a range of pulse durations, 
yielding an estimate of the minimum time required to 
implement the given operation. 
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